On the Origin of Peak-dip-hump Structure in the In-plane Optical Conductivity of the 
High Tc Cuprates; Role of Antiferromagnetic Spin Fluctuations of Short Range Order 
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An improved U(l) slave-boson approach is applied to study 
the optical conductivity of the two dimensional systems of an- 
tiferromagnetically correlated electrons over a wide range of 
hole doping and temperature. Interplay between the spin and 
charge degrees of freedom is discussed to explain the origin 
of the peak-dip-hump structure in the in-plane conductivity 
of high Tc cuprates. The role of spin fluctuations of short 
range order(spin singlet pair) is investigated. It is shown 
that the spin fluctuations of the short range order can cause 
the mid-infrared hump, by exhibiting a linear increase of the 
hump frequency with the antiferromagnetic Heisenberg cou- 
pling strength. 

PACS numbers: 74.20.Mn, 74.25.Fy, 74.25.Gz, 74.25.-q 



High Tc superconductors are the systems of strongly 
correlated electrons which show two dimensionality in 
charge transport. Various levels of gauge theoretic slave- 
boson approach to t-J Hamiltonian have been proposed 
to study high Tc superconductivity &7|. Recently we 
proposed an SU(2) slave-boson theory J6] which incorpo- 
rated coupling between the charge and spin degrees of 
freedom into the Heisenberg term. The predicted phase 
diagram showed an arch-shaped bose condensation line 
in agreement with observation. Using an improved U(l) 
slave-boson theory over our earlier one J?J , in this paper 
we study the cause of peak-dip- hump structures of ob- 
served optical conductivity ||S|-|l2fl . Various theories have 
been proposed to explain the cause of the peak-dip-hump 
structure in the optical conductivity |13|-|15||. However, 
most studies have been made to a limited range of hole 
doping and temperature, based on empirical parameters 
deduced from measurements such as the inelastic neu- 
tron scattering(INS) and the angle resolved photoemis- 
sion spectroscopy(ARPES) data. 

Using the nearly antiferromagnetic Fermi-liquid the- 
ory, Stojkovic and Pines |l3[ reported a study of normal 
state optical conductivity for optimally doped and over- 
doped systems. They showed that the highly anisotropic 
scattering rate in different regions of the Brillouin zone 
leads to an average relaxation rate of the marginal Fermi- 
liquid form. Their computed optical conductivity agreed 
well with experimantal data for the normal state of an 
optimally doped sample. Using the spin-fermion model 
jl6],[l7| and spin susceptibility parameters obtained from 
INS and NMR, Munzar, Bernhard and Cardona Q cal- 
culated the in-plane optical conductivity of optimally 



doped YBCO. Their study showed a good agreement 
with the observed peak-dip-hump structure at optimal 
doping. From the computed self energy they showed that 
the hump is originated from the hot quasiparticles and 
the Drude peak, from the cold quasiparticles. Haslinger, 
Chubukov and Abanov Jig] reported optical conductivi- 
ties <j{ui) of optimally doped cuprates in the normal state 
by allowing coupling between the spin-fermions and the 
bosonic spin fluctuations. They found that the width of 
the peak in spectral function Ak(u>) scales linearly with 
oj in both hot and cold spots in the Brillouin zone and 
cr(u>) is inversely linear in oj up to very high frequencies. 

Various studies have been limited to a restricted range 
of hole doping and temperature, relying on empirical pa- 
rameters deduced from INS and ARPES. It is thus of 
great interest to resort to a theory which depends least 
on empirical parameters and fits for a wide range both of 
hole doping(including the important region of underdop- 
ing) and temperature (encompassing the pseudogap phase 
and the superconducting phase). For this cause we use 
an improved slave-boson theory of the t-J Hamiltonian 
H that we developed recently. 

Here we briefly discuss the slave-boson theory to dis- 
cuss the coupling between the spin and charge degrees of 
freedom Q . The t-J Hamiltonian in the presence of the 
external electromagnetic field A is written, 



H=-tJ2 (e iA «cl~c ja +H.C.) 

<i,j> 

+j ( s * • s j - j n i n j) -^Yl 



<*,j> 



with Si = \^Zaf3 c ia a a0 c iP- Here is tne exter- 
nal electromagnetic vector potential; cw(cj CT ), the elec- 
tron annihilation(creation) operator at each site and a a f3, 
the Pauli spin matrix. Rewriting the electron oper- 
ator as a composite of spinon(/) and holon(&) opera- 
tors, Ci a = fio-b\ with the single occupancy constraint, 
b\bi + ^2 fjpficr — 1, we obtain the partition function, 



Z 



VfVbVXe 



- f drC 
Jo 



(2) 



with C = £\ CT //„<;./,„ • V bfdrbi + Ht-j where H t _j 
is the U(l) slave-boson representation of the above t-J 
Hamiltonian(Eq. (|J)) , 

H t -j = -t (e^flf^h + c.c.) 

<hj> 



1 



<ij> 



^ E + * E M/lJW + - 1)- (3) 



This Hamiltonian can reaily be derived from the SU(2) 
theory §. 

From the Hubbard-Stratonovich transformations in- 
volving hopping, spinon pairing and holon pairing orders 
we obtain the partition function, 



Z = J VfVbV X VA f VA b VXe So drC ' 



ss 



(4) 



with C e f f = C f + C h + Co is the Lagrangian where & = 

Si.o- fi,a®Tfi,o ' 4 ^ ^2<i,j> [xijflafj,<? + H.C.\ — 

E«j> {A( (/u/i,T " + H.C) for the 

spinon sector, C b = J2i b l d rh-tY,<i,j> {e %Aii XiMbj + 

H.C.j for the holon 
{|A{/ + i| Xl ,| 2 + 



<i,3> 



ff.cr.}-f£ <w> |Aj|»{ 

sector and £° = J(l - <5) 2 £ 

+ jE<«> |Af/|A^| 2 . Here X , Af and A fc are 
the hopping, spinon pairing and holon pairing order pa- 
rameters respectivity. 

We obtain the optical conductivity <t{uj) and the cur- 
rent response function H(lo) of an isotropic 2-D mediun in 
the external electric field E(w) by evaluating the second 
derivative of the free energy with respect to the external 
vector potential A, 

d 2 F 



<»>) = 



8J x (uj) 



dE x (u 



iu 8A X 2 



U xx (uj) 



(5) 



where J x is the induced current in the x direction; 
F = —ksThiZ ', the free energy and H xx — - "' ' 



A x =0 



the current response function in the x-direction. The 
total response function, II = II P + II D is the sum 
of the paramagnetic response function given by the 
current-current correlation function H p \r' —r,t' — t) =< 
jx(r',t')j x (r,t) > - < j x (r',t') >< j x (r,t) > with 
the current operator j x (r,t) — it(cl +xa (t)c r ^(t) — 
cj. a (t) c r + Xt a(t)) and the diamagnetic response function 
associated with the average kinetic energy, Il D =< 

K xx >= {-tT.iM^+^ + H-C-)) S 

The phase difference per unit lattice spacing associated 
with the hopping order parameter Xij = \Xij\ e<lij defines 
the gauge field, a.^ = dijO — di — 6j. The gauge fluctua- 
tions allow the back-flow condition in association with an 
interplay between the charge and spin degrees of freedom 
orriginated from the effective kinetic energy term of the 
t-J Hamiltonian. The effects of spin degrees of freedom 
are manifested through the antiferromagnetic spin fluctu- 
ations which appear in the Heisenberg exchange coupling 
term. The antiferromagnetic spin fluctuations of short 
range order(spin singlet pair) occur through the pres- 
ence of correlations between adjacent electron spins. We 



consider both the amplitude fluctuations of the spinon 
pairing(spin singlet) order parameter |A^| and the gauge 
field fluctuations. We first integrate out the spinon and 
holon fields and take the saddle point value with respect 
to the holon pairing order parameter, spinon pairing or- 
der parameter phase, the amplitude of hopping order pa- 
rameter and the Lagrangian multiplier fields in Eq.(0). 
We then obtain, 



F[A] = -fc s Tln / VfVbV x VA f VA b V\e 



-So dT ( Cf+Cb+C °) 



3 T\nJ 2?a2?|A^|(e^( F/[a ' |A/| 



+F b [A,a,\Af\[ 



„-Fb[a,|A 



(6) 



f r - f drC' 

where F 1 = -k B T\n J Vfe Jo ' is the spinon free 

energy, F b = —ksThx J Vbe So drC , the holon free en- 

rP drC 

ergy and Fq = — fc^Tln e Jo r ". The external electro- 
magnetic field couples only to the holon field but not to 
the spinon field. 

Considering the gauge and antiferromagnetic spin fluc- 
tuations up to second order we obtain the current re- 
sponse function, 



n = 



uf + u b 



n 



Ka±eL rib 

n b +nf 11 



Z TF+TTT V AA + U AA + 11 



(7) 
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where n^(II b ) is the spinon(holon) response function 
associated with the gauge field a(a and A); U XY = 

— qx^y (njyy = — Q X F dY ) i the spinon(holon) response 
function associated with both the gauge fields and the 
spinon pairing field and II AA , the response function asso- 
ciated with the spinon pairing field. It is shown that the 
first term represents the Ioffe-Larkin rule for the cur- 
rent response function contributed only from the gauge 
field fluctuations, and the second term, from the spin 
fluctuations. Each contribution comes from the coupling 
between the charge and spin degrees of freedom, as man- 
ifested by Eq.(0). 

The response function II^(n b ) is contributed from 
both the paramagnetic and diamagnetic parts. The 
paramagnetic response function is obtained from 
the current-current correlation functions Ill x (P) =< 
jl{r',t')jf{r,t) > - < j£(r',t') >< g(r,t) > for 
the spinon and ll b xx {P) =< j b (r> ' ,t')j b x (r,t) > 
j b (r',t') >< j x {r,t) > for the holon, and the diamag- 
netic response function involves the average kinetic en- 
ergy of spinon (holon). n^ A (II b A ) is given by the correla- 
tions between the spinon(holon) current and the 'anoma- 
lous' spinon(holon) pairing, II^ A =< jf(r',t')D^(r,t) > 

- < jf(r>,t>) >< Df\ r ,t) > (U b aA =< 
j b x (r',t')D b (r,t) > < j b (r>,t>) >< D»(r,t) >) 
with Df(r,t) ~ Ei{(e- iT fr,l(t)fr+i,l(t))+H.C.) and 
D b {r 7 t) - Y,i ( b r(t)b r +i(t) + H.C). Here I represents 



2 



nearest neighbor sites around location r and r = ±^(+(- 
) for x(y)-direction) is a phase to represent the spinon 
pairing of d-wave symmetry. IIaa represents correlations 
between pairing currents; n AA =< D } \r' ,t')D^ (r,t) > 
— < D* {r 1 ,t') >< Df(r 7 t) > for the spinon pairs 
and II AA =< D b (r',t')D b (r,t) > - < D b (r',t r ) >< 
D b (r,t) > for the holon pairs. 

Fig.|l| shows computed optical conductivities from the 
U(l) slave-boson t-J Hamiltonian(Eq.(||)) with J=0.3t 
for the underdoped(<5 = 0.05), optimally doped(<5 = 0.07) 
and overdoped((5 = 0.1) regions. Compared to the 
present U(l) result of optimal doping the SU(2) slave- 
boson theory || predicted a more realistic value of opti- 
mal doping close to S ps 0.15, by yielding a phase diagram 
of showing an arch shaped bose condensation tempera- 
ture in better agreement with observation. To avoid com- 
plexity, we resort to the simpler case of U(l) as our prime 
interest lies in the investigations of the role of spin fluc- 
tuations and the coupling between the charge and spin 
degrees of freedom on the formation of peak-dip-hump 
structures, since the accurate SU(2) theory will not alter 
physics on the cause of the peak-dip-hump structure. Al- 
though not shown here for other values of J we find quali- 
tative agreements with experiments in that the peak-dip- 
hump structures are well predicted below T* and Tc- In 
Fig.^| the hump peak position is seen to remain nearly 
constant with the variation of hole doping and tempera- 
ture below T* but not so above Tc- In general, the pre- 
dicted hump position tends to shift to a lower frequency 
with increasing hole concentration and with temperature, 
showing a gradual disappearance of the hump. A trend 
of rapid drop in a high frequency region is seen to be 
unrealistic. In order to find the role of spin fluctuations, 
we neglected the second term in Eq.(0). The hump struc- 
ture(dotted line in Figj^) completely disappeared, clearly 
indicating that spin-spin correlations or spin fluctuations 
associated with the spin singlet excitations are respon- 
sible for the hump formation in the optical conductiv- 
ity (Fig.||). For an additional analysis of spin fluctuation 
we computed the optical conductivity using the Lanczos 
exact diagonalization method for a two hole doped 4x4 
lattice by introducing various Heisenberg antiferromag- 
netic coupling strength J. Despite the finite size effects 
an irregular but gross feature of peak-dip-hump structure 
is still predicted indicating that the hump is originated 
from the spin-spin correlations. A linear increase in the 
hump position with J is predicted. From both the slave- 
boson and Lanczos calculations we note that the peak 
locations of the hump are sensitive to the variation of 
the antiferromagnetic coupling strength J, by showing a 
linear increase. Further, as mentioned above the neglect 
of the spin fluctuations (the second term in Eq.(Q)) led to 
a sudden disappearance of the hump structure. 

Although not shown here, using the present U(l) slave- 
boson theory the predicted spectral functions around 
(-7T, 0) point in momentum space also showed the peak- 
dip-hump structure consistent with ARPES data. This 



incoherent background or the hump around the (ir, 0) 
point was found to occur as a result of the antiferromag- 
netic spin fluctuations, having a common feature with 
the hump structure of the optical conductivity. Thus we 
conclude from these multifaceted studies that the spin- 
spin correlations or the spin fluctuations involved with 
electrons around the (n, 0) point in momentum space are 
definitely the prime cause of the hump structures below 
T* and T c . 

In the present study, by paying attention to a wide 
range of both hole doping (under-, optimal and over- 
doping) and temperature(T < Tc, Tc < T < T* , and 
T* < T) with no empirical parameters obtained from 
measurements, we examined the optical conductivity as 
a function of frequency for the two-dimensional systems 
of strongly correlated electrons. Allowing the coupling 
between the spin and charge degrees of freedom as man- 
ifested in Eq.(Q), the peak-dip- hump structures are pre- 
dicted in agreement with observations. It is shown that 
the antiferromagnetic spin fluctuations of short range as- 
sociated with the spin singlet pair excitations are impor- 
tant in yielding the observed hump structure, and that 
the hump position linearly increase with the antiferro- 
magnetic Heisenberg coupling strength. In general, the 
predicted peak-dip-hump structures are in good agree- 
ment with observations particularly in the temperature 
ranges of T < T c and T c < T < T* for the underdoped 
case. It is shown that the spin fluctuations of the short- 
est possible antiferromagnetic correlation length(that is 
the spin singlet pair) alone can cause the formation of 
the hump structure. However, considerations of both the 
antiferromagnetic spin fluctuations of correlation lengths 
larger than the spin singlet pair and the direct channel 
single spin fluctuations at high energies may be needed 
to remedy quantitative discrepancies in the rapid drops 
of optical conductivity at temperatures above T* and at 
frequencies beyond the peak location of the hump. 
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FIGURE CAPTIONS 

FIG. 1 Computed optical conductivities as a function of 
temperature for <5 = 0.05(under doped), <5 = 
0.07(optimally doped) and S — 0.1(over doped) 
cases with the antiferromagnetic Heinsenberg cou- 
pling strength of J =0.3 for all cases. 

FIG. 2 Temperature dependence of hump position as a 
function of antiferromagnetic coupling J and hole 
concentration. 

FIG. 3 The total optical conductivity (solid line) vs. a 
partial one (dotted line) contributed only from the 
first term and thus from the neglect of the spin 
fluctuation (second) term in Eq.(|7]). 
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